Modified perturbation theory for angular 
distribution in VF-pair production 



M. L. Nekrasov 

Institute for High Energy Physics, 142281 Protvino, Russia 

Abstract 

We examine the capabilities of the modified perturbation theory (MPT) for 
description of VF-pair production and decay in e + e~ annihilation. In a model 
with Dyson-resummed propagators of unstable particles, we calculate even and 
odd contributions to the distribution in the cosine of the W~ production angle 
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relative to the e~ beam. On comparing the results of calculations in the NNLO 
approximation of MPT with the exact results in the model, a coincidence of 
outcomes at the ILC energies is detected at the per-mille level. 



q : 1 Introduction 

Among the propositions for research at International Linear collider (ILC) [T], a 
significant place is assigned to measurements of interactions amongst gauge bosons. 
The problem of particular importance is the search for anomalous contributions to 
triple gauge couplings. The latter topic could become the most important issue in high 
energy physics if a light Higgs and supersymmetric particles would not be discovered. 
The primary tool for measuring the triple gauge boson vertices is the study of 
^ ! angular distributions of W bosons in the processes of VT-pair production. At the ILC a 
few million of W pairs is assumed to be produced at y/s = 500 GeV and 800 GeV [Tj[2]- 
This implies that the cross-section and various integral characteristics of angular 
distributions of W bosons may be determined at the per-mille level. (Under the 
integral characteristics we mean quantities used in the method of optimal observables 
[3111] or x 2 fit over a wide range of the phase space [5].) The anomalous contributions 
to triple gauge couplings are expected to be determined with similar precision (2J [3J 
HI E]. So, the theoretical calculation of the angular distribution must be made with 
the per-mille precision, too. This implies the next-to-next-to-leading order (NNLO) 
accuracy of the calculations, at least, on the resonant contributions of W bosons. 
Simultaneously the strict observance of the gauge cancellations must be maintained. 
The latter requirement in particular is critical in view of large gauge cancellations 
owing to which the M^-pair production is notably sensitive to triple gauge couplings. 

Unfortunately, there is a lack of calculation methods satisfying all mentioned 
requirements. For instance, the double pole approximation (DPA) successfully applied 
at LEP2 [6] can provide only the NLO accuracy of calculations and in the resonant 
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region of the cross-section only. The complex mass scheme (CMS) [7J provides the 
NLO accuracy at higher energies, but the NLO is not sufficient for support of all 
potentialities at ILC. So a more powerful scheme is required which would allow one 
to consistently involve the NNLO in the calculations. 

A probable candidate for such scheme is a modified perturbation theory (MPT) 
[8j HI [10] • Its main feature is the direct expansion of the probability instead of ampli- 
tude in powers of the coupling constant with the aid of distribution-theory methods. 
The latter methods permit to impart a well-defined meaning to the resonant contribu- 
tions of unstable particles in the expanded probability. A condition of asymptoticity 
(and therefore of completeness) of the expansion must ensure the gauge cancellations. 
In the case of pair production of unstable particles, the most-elaborated description 
of the method is given in [10]. In [UJ [12] the convergence properties of the MPT ex- 
pansion for the total cross-section were tested in models related to the top-quark-pair 
and W-pair production. (The models were based on the "improved" Born approxi- 
mation and the Dyson resummation up to three loops of the resonant contributions of 
unstable particles.) In both cases a good convergence of the MPT series was detected 
at the ILC energies. In the case of the TV-pair production, the precision of the NNLO 
approximation was observed at the level better than one per-mille |12j . 

In this paper we consider angular distribution of the W~ boson relative to the 
direction of the e~ beam in the center-of-mass frame. Actually this distribution is used 
in the experimental analysis of the triple gauge couplingslj In fact the distribution 
is analysed together with changing the beam polarization in order to disentangle 
contributions from ^WW and ZWW vertices. From standpoint of the MPT, the 
changing of the polarization manifests itself as the changing of test function. However 
a variation of the test function has virtually no effect on the convergence of the MPT 
series [TJJ, [12] provided that the set of singularities arising in the MPT is unchanged. 
The set of singularities is invariable in the cases of even and odd contributions to 
the distribution in the cosine of the W production angle. So for the analysis of the 
convergence properties of the MPT series, we may confine ourselves to consideration 
separately of the even and odd contributions and in the case of unpolarized beams 
only. By these means we would be able to verify the convergence properties of the 
MPT expansion simultaneously for wide variety of beam polarizations. 

The paper is organized as follows. In Sec. [2j we consider specifics of the MPT in 
the case of the angular distribution and determine a method of handling new-type 
singularities that arise in this case in the framework of the MPT. In Sec. [31 we present 
results of numerical calculations. In Sec. [U we discuss results and make conclusions. 

2 MPT for angular distribution 

The angular distribution of W^-pair production in e + e~ annihilation has the form of a 
convolution of hard-scattering contributions with the flux function describing initial 

^dn the case of transverse beam polarization also the azimuthal-angle distribution of W bosons 
is considered [HE], but it does not contain new MPT entities. For this reason we do not discuss it. 
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state radiation, 



r ds' 

Va(s,cos9) = — <j)( s' /s;s)Va(s', cos 9) . (1) 
J* ■ s 

Here T> = d/dcos#, 9 is the W~ production angle relative to the e~ beam in the rest 
frame, <fi is the flux function. The hard-scattering contributions T>a are described as 
an integral of exclusive angular distribution over the virtualities of W bosons, 

poo poo 

Va(s',cos9) = / ds 1 ds 2 (1 + 5 soft ) Va ex (s\ cos 9 ; s u s 2 ) . (2) 

Slmin J *2min 

In this formula Si min and s 2m in are the minimum virtualities, ^Js~^ l = ^Si min + 
v /s 2m in, 5 so ft stands for factorized contributions of soft massless particles. The exclu- 
sive contribution, in turn, is written as a product of Breit-Wigner (BW) factors p(sj), 
some kinematic factors, and a function $ which is the rest of the amplitude squared, 

& ex (s,cos9;s 1 ,s 2 ) = &(y/s-y/sl-y/si)y/\(s,s 1 ,s 2 ) $(s,cos0; si,s 2 ) p(s 1 )p(s 2 ) ■ (3) 

Here X(s,si,s 2 ) = (s —Si — s 2 ) 2 — 4siS 2 , and 0(. . .) is the step function. Function 
$ corresponds to one-particle irreducible contributions in the channels of W bosons. 
Therefore it has no singularities on the mass-shell. As a first approximation, we 
regard that $ is an analytic function and has no singularities. (By this we mean that 
the appropriate singularities are to be considered on the fact they are identified.) On 
the contrary, the kinematic factor, which is the product of the step function and the 
square root of the kinematic function, explicitly has singularities. The BW factors, if 
we naively expand them in powers of the coupling constant, generate non-integrable 
singularities. 

However, the singularities might be made integrable if we expand the BW factors 
in the sense of distributions [TB]. Actually, this is a basic idea of the MPT approach 
[8] . In this case the expansion of a separately taken BW factor is beginning with the 
5-function, which corresponds to the narrow-width approximation. The contributions 
of the naive Taylor expansion are supplied with the principal-value prescription for 
the poles. The nontrivial contributions are the delta-function and its derivatives with 
coefficients C n , which are polynomials in the coupling constant a and which correct 
the contributions of the Taylor. Within the NNLO, the expansion looks as follows [TO] : 

p(s) = ^ \s- A^ 2 + S(s)r 2 = 5(s-M 2 ) (4) 

2 

+ Y,Cn(a)5^(s-M 2 ) + 0(a 3 ). 

n = 

Here PV is the principal-value prescription, M is the renormalized mass of the un- 
stable particle, T is its Born width, aSi(s) is the one-loop self-energy. Within 
the NNLO in OMS-like schemes of UV renormalization, coefficients C n (a) include 
three-loop self-energy contributions and their first derivatives determined on-shell. 
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Unfortunately, the above expansion is senseless if the weight at the expansion is not 
smooth enough. In formula ([3]) the weight is not smooth because of singular behavior 
of the kinematic factor. So, the next ingredient of the MPT is the analytic regu- 
larization of the kinematic factor via the substitution [A(s, s\, s 2 )} 1 ^ 2 — > [A(s, s\, s 2 )]^ 
|10j . The rest of the weight is represented in the form of Taylor expansion in variables 
Si and S2 around the mass-shell, with a remainder. The integrals of the remainder 
may be numerically calculated without the regularization. The integrals of Taylor are 
reduced to the sum of basic integrals, which may be analytically calculated irrespec- 
tive of details of definition of the weight. In the end with v — 1/2 this gives finite 
outcomes in any order in the expansion in powers of a and the expansion remains 
asymptotic [10J. 

When calculating the total cross-section with $ determined in the Born approxi- 
mation, this scheme works well [12]. However in the case of the angular distribution 
a non-analyticity in $, related to neutrino propagator in the t-channel, becomes rele- 
vant. Namely, the square root of the kinematic function A in the denominator of the 
propagator becomes relevant, 

(5) 



s — si — S2 — y/\ cos 9 t u 

Note that in the total cross- sect ion, the integrating dcos 9 results in an entire function 
of A in the vicinity of A = 00 But in the case of angular distribution the a/A persists. 
Of course, the y/X may be analytically regularized by means of the substitution \/\ — > 
\ u . However, when calculating the Taylor series, this leads to a growth of singularity, 
which considerably complicates calculations. 

An alternative way is to represent propagator ([5]) in an equivalent form with an 
entire function in the denominator, 

s - si - s 2 + y/X cosfl 
" ~ (s- Sl -s 2 ) 2 -X cos 2 9' 1 ' 

Then, function $ may be considered as a sum of two contributions, one with integer 
powers of A and another with an additional factor a/A in the numerator. The former 
contribution is handled by the above scheme. In the latter contribution we remove 
the additional y/~X in $ in favour of the kinematic factor, so that the kinematic factor 
would include \ v+1 / 2 instead of \ u . 

Further let us remember that the contributions of the above-mentioned Taylor 
expansion are expressed in terms of the basic integrals 1%, J 1 ^-, ^ lBs , 5^ i ri2 , and 
C" n2 . (In fact C" n2 is represented as a sum of other basic integrals and a regular 
function. For this reason we do not discuss C% n2 below.) In the original work [10] the 
basic integrals were calculated at non-integer u, and only half-integer v was needed 
in the final stage with smooth $. At the same time, a transition to integer v was 
not obvious because of singular behavior of some factors. However, the transition is 
possible through asymptotic expansion of singular factors. The mentioned factors are 



2 Although the explicit result at first glance includes singularity at A = 0, see e.g. [14 , it is 
represented as a Taylor in A in the vicinity of A = 0. 
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the T-f unction whose argument tends to negative integer or zero, and the distribution 
x+ with f3 tending to negative integer. Their asymptotic expansions are as follows: 

r(i - jv+£)= ^(S +o(i) ' (7) 



-N+e 



1 f- 



,JV-1 



SW-Vfx) + xl N + 0(e) 



e (JV-1)! 

Here N = 1,2, ... , and x+ N is the adjoint distribution of the 1-st order [13]. Fortu- 
nately the x+ N does not appear in the final formulas owing to the relation 

x~ N + (-y N (-x)~ N = PVx~ N . (9) 

It is worth to recall that the asymptotic expansion in (J5]) is defined in the sense of 
distributions, which means that both sides of the relation are to be integrated with 
smooth enough weight. 

Formulas ©-(E]) are sufficient for the determination of the basic integrals with 
non-negative integer v. The basic integrals needed within the NNLO are as follows: 

\\{x)=x + , (10) 

I ^ I 

J 1 (x\ai) = — (x + di) + + Q(x + di) x In — , (11) 

a, 

1 \x\ 
J 2 (x\ai) = (x + cti) + — Q(x + di) In — , 

Al,=x + , Al 2 = A^ = -Q(x), Al z = A 3 \ = i 8{x) , A 2 \ = 5(x) , (12) 
B^xlai) = Bl^xldi) = Jl(x\ai), (13) 



B 2 \(x\ai) = 5 2 \(a;|ai) = J 2 (x\ai) . 

Here x is dimensionless variable associated with s by the relation x = 2M~ 1 ( v /5— 2M), 
and dimensionless parameters, = 2M 1 (M — y/si^). 

Another problem arising at calculating the angular distribution is connected with 
zeros in the denominator in Actually, the denominator is nulled at s — Si — s 2 = 
and cos 6 = 0. The former equality implies A < 0, which means non-physical region. 
However, the Taylor expansion of the weight effectively means extension into non- 
physical region. Really, with small enough s one of Sj, or both s± and S2, are taken 
due to the Taylor on the mass-shell. As a result, the denominator in ([6]) affects the 
calculations with s — si — s 2 = 0. This means failure of the calculations at cos 6 = 
and an instability of outcomes of calculations in the vicinity of cos 8 = 0. 

The problem may be solved by redefinition of the denominator in the non-physical 
region so that it would be positive everywhere. So, by taking advantage of the identity 

(s-si- s 2 f - X cos 2 9 = A (1 - cos 2 9) + A Sl s 2 , (14) 
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we represent fl6]) in the form 

s - si - s 2 + y/\ cos 6* 
" ~ #(A)(l-cos 2 fl)+4 Sl s 2 1 j 

with 

H(\) = 0(A) A + 9(-A)/i(A). (16) 

Here h(X) is a function that provides positivity of the denominator and, simultane- 
ously, sufficient smoothness of the function H(X). Since within the NNLO the Taylor 
implies second-order derivatives, function H(\) must be continuous with its second 
derivatives. These conditions are satisfied, for instance, by 

h(\) = A [1 - (-A/Ao) 71 ] , (17) 

where n > 2, < A < 4 si min s 2min (n+l)( n+1 )/ n /n. 

We stress that the above procedure in no way affects the integrand ([2]) in the 
region of integration (in the physical region) and thus in no way affects the value 
of the integral itself. Nevertheless, the redefinition of the integrand in the extended 
region allows us to calculate the asymptotic expansion of the integral by means of 
asymptotic expansion of the BW factors and the Taylor expansion of the weight. The 
result of the calculation must be the same as that in the case of direct asymptotic 
expansion of the previously calculated integral. 

Now we briefly discuss the physical model in which we carry out calculations. As 
mentioned in the Introduction, the model is the same that was used in [12]. It is based 
on the improved Born approximation in the Standard Model for the amplitude and 
the Dyson resummation up to three loops of the resonant contributions of unstable 
particles. The improved Born approximation implies taking into consideration the 
universal soft massless-particles contributions that are collected in the flux function 
and in the Coulomb factor in 5 so ft. We consider the flux function in the leading- 
log approximation, and the Coulomb factor in the one-photon approximation with 
specific resummation that does not affect the BW factors, see details and references 

in dm eng. 

The scalar part of the IV-boson propagator in the model, we consider in the form: 

A _1 (s) = s — M 2 + aReSi(s) + i«ImSi(s) 

+ a 2 [R 2 + i/ 2 + (s-M 2 )(R' 2 + H' 2 )] + a 3 (R 3 + H 3 ). (18) 

Here R n = Re£ n (M 2 ), I n = Im£ n (M 2 ), and the prime means derivatives. ReEx(s) 
and ImS^s) can be determined by direct calculations, but for simplicity we restrict 
ourselves by considering quark and lepton contributions only. In this way we avoid 
the IR divergences. The on-shell self-energy contributions can be determined by 
means of the UV renormalization conditions and the conditions of unitarity. We 
use the OMS scheme [T5J [16] of the UV renormalization, a natural generalization 
of the conventional OMS scheme for the case of unstable particles. In this scheme 
the imaginary contributions to on-shell self-energies are expressed through the Born-, 
one-loop, and two-loop contributions to the widths of the unstable particles. 
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Propagator (ITS]) is considered as the exact one in the model. The MPT-expansion 
of the appropriate BW factors are determined by (J4j) with coefficients C n (a) given 
in [TU]. The definition of $ in the Born approximation is straightforward except Z- 
boson propagator in the s-channel. Namely, at the resonant energies for the Z boson, 
we define its propagator as was done in f JTSj) . So, when considering the model in the 
MPT-mode, we expand Z-boson propagator in the MPT sense. However, at high 
enough energies, we define Z-boson propagator to be free, thus preserving the gauge 
cancellations and the unitarity, see details in [12]. 

The angular distributions in the model are straightforwardly calculated. We calcu- 
late the case of unpolarized beams. The outcomes with the propagator (fT8|) . are called 
the "exact results in the model". The angular distributions in the MPT-mode have 
the form Va = Va + Vaai+Va 2 a2, where o"o is the cross-section in the LO approxi- 
mation, aa\ and a 2 a 2 are the NLO and NNLO corrections. So, T><r NLO = Va +Vaai 
and Vcjnnlo — Va + Vaai + Va 2 a 2 determine the NLO and NNLO approximations, 
respectively. 

3 Numerical results 

For numerical calculations, we use the same values of kinematic and model parameters 
and the same parametrization that was used in [12] . The outcomes of the calculations 
are presented at the typical ILC energies y/s = 500 GeV and yfs = 800 GeV [TJ [2] . 

In Figs. [THE] we present the full cos ^-distribution, and separately the even and 
odd contributions, respectively. In all figures the thick curves show the exact result 
in the model. The dotted, dashed, and thin continuous curves show the results in the 
LO, NLO, and NNLO approximations in the MPT, respectively. The curves for the 
NLO and NNLO are omitted in the upper row in Fig. [TJ since they merge with the 
curves for the exact result at a given scale. Instead, they are shown in the middle 
and lower rows, where the percentages with respect to the exact result are presented. 
In Fig. |2] and El we omit "upper row" with the absolute- value results for the even and 
odd contributions in view of their minor informativity, and present only percentages 
of the corresponding approximations relative to the appropriate exact results. In all 
figures the lowest row repeats the results of the preceding row with greater scale on 
vertical axis. 

Examining the figures, we note first and foremost a stable behavior of the NLO ap- 
proximation and more stable behavior of the NNLO approximation in the cases of even 
and odd contributions (Figs. [2]and[3j respectively). The accuracy of the NLO approx- 
imation in these cases is within a few per-mille and that of the NNLO approximation 
persistently is better than one per-mille. In the case of the full angular distribution 
(Fig. [lj the picture generally is the same except the region of cos 9 approaching —1. In 
the latter region the even and odd contributions almost cancel each other, thus giving 
in the sum a substantially reduced outcome. For instance, with cos# = —0.95 the out- 
come is diminished by about two orders of magnitude: Va = 12.72 — 12.59 ~ 0.12 pb 
and Va = 5.574 - 5.513 = 0.061 pb for ^Ts = 500 GeV and 800 GeV, respectively. 
As a result, a discrepancy of the sum increases near cos 6 = —1. The tendency is 
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Figure 1: Distribution in the cosine of the W production angle relative to the e 
beam in e+e" ->• 7, i? — »■ W + W~ 4/ at a/s = 500 GeV (l.h.s.) and y/s = 800 GeV 
(r.h.s) in the case of unpolarized beams. In the upper row, the exact result in the 
model and the result in the LO in the MPT are shown. At the middle and lower 
row, the percentages of the results in the LO, NLO, and NNLO approximations in 
the MPT with respect to the exact result are shown. 
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Figure 2: The MPT results in percentage with respect to the exact result for the even 
contribution to the cos ^-distribution. 
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Figure 3: The MPT results in percentage with respect to the exact result for the odd 
contribution to the cos ^-distribution. 
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Table 1: Results of MPT calculations for forward-backward asymmetry Afb- 



y/s (MeV) 


A EXACT 
A FB 


a lo 

FB 


aNLO 
FB 


aNNLO 
FB 


500 


0.9025 


0.9150 


0.9011 


0.9022(1) 




100% 


101.38% 


99.94% 


99.96(1)% 


800 


0.9207 


0.9409 


0.9202 


0.9204(1) 




100% 


102.19% 


99.95% 


99.97(1)% 



as follows: with cos 6* = 0.95 at y/s = 500(800) GeV the discrepancy in the NNLO 
approximation makes up 0.05(0.03)%; with cos 9 = — 0.64(— 0.63) the discrepancy 
overcomes the level of 0.1%; with cos 6* = —0.95 the discrepancy makes up 1.4(1.5)%. 
The latter values are to be compared with the discrepancies for the even and odd 
contributions taken separately, which are 0.04(0.02)% and 0.05(0.04)%, respectively, 
and with two orders of magnitude cancelled in the full cos ^-distribution. 

It should be remembered that with cos 9 approaching —1 the quality of the experi- 
mental data in the case of unpolarized beams becomes worse, too. This occurs because 
of the smallness the values of the angular distribution and therefore its statistical sig- 
nificance decreases, and also because of the increasing background [5]. So the effect 
of decreasing the accuracy of the MPT description most likely will not spoil quality 
requirements in this region. On the other hand, with polarized beams the angular 
distribution must become more noticeable at cos 9 ~ — 1. In particular, with purely 
right-polarized e~ beam the odd contribution disappears. Therefore the very reason 
for the above mentioned cancellations disappears, too. So, in general, in most cases 
of beams polarization there should be no peaks in the region of backward-scattering 
in the percentage distributions of the MPT description. 

In the end of this Section it is worth noticing that the above mentioned effect of 
decreasing accuracy of the MPT description in the case of unpolarized beams and 
large backward scattering angles does not appear at all in the cases of the total cross- 
section and the forward-backward asymmetry Afb- Really, the total cross-section 
is determined as the integral dcos^ of the cos ^-distribution in symmetric limits, so 
only the even contribution is relevant. Therefore, the accuracy of the description of 
the total cross-section in the NNLO is better than one per-mille at the ILC energies 
[j~2] . The forward-backward asymmetry A ps is determined as the ratio of the integral 
of the odd contribution to the integral of the even contribution to cos ^-distribution. 
Therefore, the accuracy of the description of Aps is within one per-mille, as well. 
The numerical results at the characteristic ILC energies are presented in Table [TJ 
The numbers in parenthesis in the last column show uncertainties in the last digits 
due to computations. In other columns the uncertainties are omitted as they are 
beyond the precision of the presentation of data. The algorithm for estimating the 
uncertainties is described in [T7] . 
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4 Discussion and conclusions 



We have tested the applicability of the MPT in the case of angular distribution of W 
bosons in the processes of VT-pair production in e + e~ annihilation. Specifically, in 
a model that admits exact solution, we have calculated separately the even and odd 
contributions to the distribution in the cosine of the W~ production angle relative to 
the e~ beam. At the ILC energies a coincidence of the MPT outcomes with the exact 
result is detected within a few per-mille in the NLO approximation, and within one 
per-mille in the NNLO approximation. 

In fact, despite the results of previous investigations of the total cross-section [12], 
the accuracy of the MPT-description in the case of angular distribution a priori was 
not known because of the presence of new-type singularities in the odd contribution to 
the cos ^-distribution. Such singularities were missed in previous MPT investigations. 
Initially they arise from the t-channel neutrino exchange and then they give rise to 
logarithmic and delta-function singularities in the exclusive cross-section, in addition 
to the power-like singularities in the general case [10j. (Recall that in the observable 
cross-section all singularities are integrated out.) Another novelty in the present 
calculations is a trick permitting to avoid spurious singularities arising in the non- 
physical momentum region involved in the MPT calculations. The mentioned trick 
enriches the set of techniques in the approach of MPT. 

Actually, we have carried out calculations in the case of unpolarized beams only. 
However, the results may easily be generalized to cases with different polarizations. 
This is possible due to intrinsic property of the MPT, which is the almost indepen- 
dence of the results expressed in the relative units from smooth variation of the test 
function (i.e. without the occurrence of new singularities). This property is quite 
natural from the standpoint of theory of distributions where the distributions are 
considered as continuous linear functional acting on linear space of test functions 
[T5] . In the MPT this property was explicitly verified [TTj [T2"] . It permits to extend 
the result concerning the accuracy of the MPT description to cases with different 
polarizations. Really, a change of the polarization means smooth variations of the 
test functions for the even and odd contributions to the cos ^-distribution. Owing to 
smoothness of the variations, the accuracy of the MPT description of the even and 
odd contributions must remain at the same level. 

An important corollary of this result is that the accuracy of the description of the 
full angular distribution must remain at the same level, i.e. at the level of one per-mille 
in the NNLO approximation. (Let us recall that this is the accuracy that is required 
at the ILC.) The only exception is the case of large cancellations between the even and 
odd contributions, what happens in the case of unpolarized beams in the region of large 
backward scattering angles. In this region the accuracy in relative units is decreasing 
following the scale of the cancellation in the full angular distribution. However, owing 
to diminishing the signal the quality of experimental data becomes worse, too. So it is 
unlikely that the mentioned effect would be crucial when extracting information from 
experimental data. At the same time, in most cases of beams polarization when there 
are no large cancellations and therefore the angular distribution is not too small, the 
accuracy of the MPT description in the NNLO must remain at the per-mille level 
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throughout. 

To conclude the discussion, two more comments are in order. First, owing to the 
almost independence of the results of MPT calculations expressed in relative units 
from particular form of the test function, the turning on of the loop correction to the 
amplitude should not lead to noticeable modification of our results concerning the 
precision of MPT calculations. Notice this remark should be in force also because 
all types of basic singularities that can appear in the hard-scattering cross-section 
have been taken into consideration in our analysis (the power-like, logarithmic, and 
delta-function). The second remark concerns the single- resonant and interference 
contributions. Actually this issue was discussed earlier [10], and it was shown that 
such contributions may be taken into consideration precisely by the same technique 
as discussed above. So in these cases our results about the precision of the MPT- 
description must be in force, as well. 

In summary, the MPT in the NNLO approximation gives on the whole satisfactory 
(from standpoint of ILC) results for the angular distribution of W bosons in the 
processes of P^-pair production and decay in e + e~ annihilation. On this basis and 
taking into consideration the preceding results for the total cross-section [12J, we 
conclude that the MPT is really a good candidate for theoretical support at the ILC 
of the processes of W^-pair productions. 
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